Exploiting quantum parallelism of entanglement for a complete experimental 
quantum characterization of a single qubit device 
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We present the first full experimental quantum tomographic characterization of a single-qubit 
device achieved with a single entangled input state. The entangled input state plays the role of all 
possible input states in quantum parallel on the tested device. The method can be trivially extended 
to any n-qubits device by just replicating the whole experimental setup n times. 
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The new field of Quantum Information Q has recently 
opened the way to realize radically new processing de- 
vices, with the possibility of tremendous speedups of 
complex computational tasks, and of cryptographic com- 
munications guaranteed by the laws of physics. Among 
the many problems posed by the new information tech- 
nology there is the need of making a complete experimen- 
tal characterization of the functioning of the new quan- 
tum devices. As shown recently in Ref. ||, quantum 
mechanics provides us with the perfect tool to achieve 
the task easily and efficiently: this is the so called quan- 
tum entanglement, the basis of the quantum parallelism 
of future computers. In this letter we present the first 
full experimental quantum characterization of a single- 
qubit device, based on this method. Since the method 
can be easily extended to any n-qubit device, the present 
experiment represents a first test of the feasibility and of 
the experimental limits of the new general tomographic 
method. 

How we characterize the operation of a device? In 
quantum mechanics the evolution of the state is com- 
pletely described by the so called quantum operation^ of 
the device, that here we will denote by E. More precisely, 
the output state p ou t is given by the quantum operation 
E applied to the input state pi n as follows 



Pout 



E (pin) 

Tr[E (p m )Y 



(1) 



The normalization constant Tr [E (pi n )] in Eq. ([!]) is also 
the probability of occurrence of the transformation E (e. 



* Electronic address 
i Electronic address 

* Electronic address 



|fraiiccsco.dcmartim(Q umromal .it 
kndrea.mazzci@uniromal.it| 



arco.ricci@uiiiroma l.it 
URL: frittp : //atesi^TTphys 



.umroma 



■ Electronic address 
**URL: 



sriann(5j)unipv,it 

tittp : / /www . qubit . itl 



it 



t^Also at Department of Electrical and Computer Engineering, 
Northwestern University, Evanston, IL 60208 



g. when there are other possible alternatives, such as 
when we consider the state transformation due to a mea- 
suring device for a given outcome). Therefore, apart from 
a normalization factor, the quantum evolution is always 
linear, with the quantum operation playing the role of 
the so-called transfer matrix of the device, a mathemati- 
cal tool very popular in optics and electrical engineering. 

Now the problem is: how to reconstruct the form of 
E experimentally? Since E is essentially a transfer ma- 
trix for a linear system, one would be tempted to adopt 
the conventional method of running an orthogonal ba- 
sis \n) of inputs and measuring the corresponding out- 
puts by quantum tomography^. However, since states 
are actually operators — not vectors — as a consequence of 
the polarization identity in order to get all off-diagonal 
(complex) matrix elements of the state one actually needs 
to run not the basis itself, but all the linear combina- 
tions of its vectors 2~2 (|n') + n\n")), with n — ±l,±i. 
In the following we will call such sets of states faithful, 
since they are sufficient to make a complete tomography 
of a quantum operation. This method is essentially the 
quantum process tomography given in Ref. 0, which was 
experimentally demonstrated in NMR § and recently in 
quantum optics . The main problem with such method 
is the fact that in most practical situations faithful sets 
of input states are not feasible in the lab. For example, 
for continuous-variables process tomography in the Fock 
representation, the states \n'} and \n") would be photon 
number states, and achieving their superpositions will re- 
main a dream for experimentalists for many years. But 
here the quantum parallelism of entanglement comes to 
help us, with a single input entangled state that is equiv- 
alent to running all possible input states in parallel [Q. 
Thus, we don't need to prepare a complete set of states, 
but just a single entangled one, a state commonly avail- 
able in any modern quantum optical laboratory. 

In the following we will use the double-ket notation 
|^)) G Tl <8> H. to denote bipartite states corresponding to 
the matrix ^ nm of coefficients on fixed given orthonormal 
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basis 



Iran) of H®H 



nm). 



(2) 



In our optical implementation the entangled systems con- 
sist of two single-mode optical beams, and the Hilbcrt 
space is two-dimensional, since we will consider only 
single-photon polarization states. The key feature of the 
method implies that only one of the two systems is input 
into the unknown transformation E, whereas the other is 
left untouched, as in Fig. [jj This setup leads to the out- 
put state Rout, which in tensor notation writes as follows 



fl^ = E® Ifltf)) ((*!), 



(3) 



where I denotes the identical operation. It is a result of 
linear algebra that R ou t is in one-to-one correspondence 
with the quantum operation E, as long as the state |\&)) 
is full-rank, i. e. its matrix * is invertible. This is the 
case of a so-called maximally entangled state, where the 
matrix \P is proportional to a unitary one. Full-rank en- 
tangled states can be easily generated by spontaneous 
parametric down conversion of vacuum, as in the exper- 
iment reported here. Note that even when a faithful set 
of input states is available in the lab — which is actually 
true in our case of single-photon polarization states — 
nevertheless a single faithful entangled state can be much 
more efficient and more practical. As a matter of fact, in 
practice, generation of single-photon polarization states 
relays anyway on entanglement, and, as we will see in 
the following, the present method uses all experimental 
data much more efficiently than the conventional quan- 
tum process tomography jl], ^| . 
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FIG. 1: General experimental scheme of the method for the 
tomographic estimation of the quantum operation E of a sin- 
gle qubit device. Two identical quantum systems — two op- 
tical beams in the present experiment — are prepared in an 
entangled state 1^)). One of the two systems undergoes the 
quantum operation E, whereas the other is left untouched. At 
the output one makes a quantum tomographic estimation, by 
measuring jointly two observables (each for each beam) from 
a quorum {Q(l)}. In the present experiment the quorum is 
represented by the set of Pauli operators a x , a y and a z . 

Now, how to characterize the entangled state R «t at 
the output? For this purpose a technique for the full 
determination of the quantum state has been introduced 
and developed since 1994. The method named Quantum 
Tomography Q has been initially introduced for the state 
of a single-mode of radiation, the so called Homodyne 



Tomography, and thereafter it has been generalized to 
any quantum system. The basis of the method is the 
measurement of a suitably complete set of observables 
called quorum. In our case we need to measure jointly 
two quora of observables on the entangled qubit, here the 
quorum being the Pauli matrices o~ x , o~ y , cr z . The qubit 
is encoded on polarization of single photons, as follows 



|Q) = /i f |0), 



ID 



»t|n), 



(4) 



where |fi) denotes the e. m. vacuum, and h, h' and v, 
v' the annihilation and creation operators of the hori- 
zontally and vertically polarized modes associated to a 
fixed wave- vector k, respectively. In synthesis, Eq. (^) 
means that the " logical zero" is encoded on a single hor- 
izontally polarized photon, whereas the "logical one" is 
encoded on a vertically polarized photon. In the present 
representation, the Pauli matrices write as follows 

o~x = h'v+v' h, o~ y — i (h'v — v'h) , a z = h^h—v^v. (5) 



According to Eq. (F 
achieved as in Fig. ] 



the <7 2 -photo-detector is simply 
In order to understand how to 
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FIG. 2: Pauli-matrix a z measurement apparatus for photon- 
polarization qubits inserted at the end of each optical beam. 
The beam is split by a polarizing beam splitters (PBS) into 
its horizontal and vertical components, which are separately 
detected and recorded with plus and minus sign, respectively. 
For measuring the other two Pauli matrices a x and a v , the 
PBS is preceded by a suitably oriented 4 and 
respectively (see text). 



j wave-plate, 



design detectors for o~ y and cr z we still need some simple 
algebra for wave-plates. The ring of Pauli matrices is 
completed by including the "identity" CTo = h'h + v'v, 
corresponding to single-photons states. In the following 
we use the popular relativistic conventions, denoting by a 
the three- vector of operators a — (p\ , 02, 03) and by a the 
tetra- vector a = (00, o~i, o~%, 03), and use Greek indices for 
three-vectors components a — 1, 2, 3 (or a = x,y,z), and 
Latin indices for tetra- vector components i = 0, 1, 2, 3. 

A wave-plate changes the two radiation modes accord- 
ing to the matrix transformation 
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where the matrix 



is given by 
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CZ- 
SZ- 



SZ- 
z+ — CZ- 



(6) 
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where s = sin 29, c = cos 29, z± = | (l ± e 1 ^), 9 is the 
wave-plate orientation angle around the wave-vector k, 
(ft = 2ir5/\, A is the wave-length, and 5 is the length 
of the optical path through the plate. Special cases are 
the j wave-plate which can be used with 9 — 7r/4 to 
give the right c+ and left c_ circularly polarized modes 
c± = e 47r / 4 2~2 [±h + iv), and the -| wave- plate which can 
be used to give the diagonal linearly polarized modes 
d± = 2~z(h ± v). On the Pauli operator vector a the 
mode transformation due to a wave-plate writes 



with rotation matrix 



S>,9 



s 2 +c 2 cos<p — c cos (ft sc(l 

c sin (ft cos (ft 
sc( 1 — cos (ft) ssmcft c 2 4 



— cos$ 
s sin (ft 
■ s 2 cos < 



(8) 



(9) 



From Eq. (|) we can see that a <7 T -detector can be ob- 
tained by preceding the ovdetector with a -| wave-plate 
oriented at 9 = 7r/8, whereas a o^-detector is obtained by 
preceding the cr z -detector with a h wave-plate oriented 
at 9 = 7r/4. When collecting data at a <r a -detector, we 
will denote by s a — ±1 the corresponding random out- 
come, s a = +1 corresponding to the /i-detector flashing, 
and s a = — 1 to the w-detector flashing. The general ex- 
perimental setup is then given by two Pauli detectors — 
for measuring o~ a and 073 for varying a and (5 — at the 
output of the entangled beams, as in Fig. [j]. The experi- 
mental data are collected in coincidence, with two of the 
four photo-detectors firing, one for each Pauli detector, 
thus guaranteeing that the result will be essentially unaf- 
fected by quantum efficiency The experimental correla- 



tions s^P sf^ of the random outcomes s) n> of the detector 

* J 1 

at the n-th beam (n = 1,2) on the entangled state I*)) 
must coincide with the following theoretical expectations 
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of } )|*)) = Tr [V+aiVa*] , (10) 



,(2) 



Sg 1 - 1 s ■ 2 ' ) . For 



and, obviously, s\ = s ■ Sq and s 
maximally entangled states we have the isotropy condi 

tion Sa — Sa = for a = x,y, z. The four Bell states 
will correspond to the four Pauli matrices o~j via a state 
coefficients matrix ^ given by ^ — -j^Vj- On the other 
hand, when a quantum device performing the unitary 
transformation U is inserted in one of the two entangled 
beams as in Fig. [j], the entangled state I*)) is changed to 
?7® /I*)), which corresponds to the new coefficients ma- 
trix * — > U^S . In our lab we used the "triplet" state cor- 
responding to * = <j x / \[2, which is generated via sponta- 
neous parametric downconversion by an optical paramet- 
ric amplifier physically consisting of a nonlinear BBO (fi- 
barium-borate) crystal plate, 2 mm thick, cut for Type II 
phase matching and excited by a pulsed mode-locked ul- 
traviolet laser UV having pulse duration r = 140 fsec and 



wavelength A p =397.5 run. The wavelength of the emit- 
ted photons is A = 795nm. The measurement appara- 
tus consisted of two equal polarizing beam splitters with 
output modes coupled to four equal Si-avalanche photo- 
detectors SPCM-AQR14 with quantum efficiencies QE 
~ 0.42. The beams exciting the detectors are filtered by 
equal interference filters within a bandwidth AA = 6nm. 
The detector outputs are finally analyzed by a computer. 

With the above apparatus we want now to experimen- 
tally determine the matrix elements of the state |*)) in 
Eq. (0). From the trivial identity 



(ran = 



(11) 



we obtain the matrix ^ nm for the input states in terms 
of the following ensemble averages 



((#|01)(nm|W)) 



(12) 



where exp (icp) — ^oi/l^oil is an unmeasurable phase 
factor. The choice of the vector |01) is arbitrary, and it is 
needed only for the sake of normalization, e. g. we could 
have used |10) or |11), instead. Using the tomographic 
expansion over the four Pauli matrices Q we see that, 
via Eq. (|l2"|), the matrix element of the input state is 
obtained from the following experimental averages 
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(13) 



where 



p = «*|01><01|*» = -(1 + 4 X) )(1 - 4 2) ) (14) 

is the fraction of events with one (T z -detector firing on h 
and the other on v, and the matrix Q\\ m is given by 



5 nm 



(n\ai\0)(m\a 3 \l} 



(15) 



The unitary matrix U nm of the device is now obtained 
with the same averaging as above, but now for the state 
at the output of the device 1*7*)) = (17®/) |*)), namely 



((*7*|01)(tito|**7)) 



(16) 



where we use again Eqs. (|14|) and (|15|), but now the 
average expressed by Eq. (jlj) is carried out over the 
output state | U *)) . The (complex) matrix elements U nm 
are obtained from Eq. (|16|) by matrix inversion. This is 
of course possible since the matrix * is invertible, due to 
the maximally entangled character of |*)). 

An experimental demonstration of the tomographic 
method is given in Figures |^ and ^, where both the real 
and imaginary parts of the four measured matrix ele- 
ments of the unitary operator U of the analyzed device 
are reported for two different devices, and compared with 
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U[* = 0.45m 6- -0.13911] 



real part 




U()() Uoi Uio U|i 



FIG. 3: Experimental characterization of a single optical 
wave-plate with retardation phase <j> = 0.45-7T and orientation 
angle of the optical axis respect to the laboratory horizontal 
direction 9 — — 0.1387T. The experimental matrix elements 
Unm of the wave-plate are superimposed to statistical errors 
for 8000 events, and compared with the theoretical values. 

the theoretical values. As one can see, the experimental 
results are in very good agreement with theory, within 
experimental errors. As a first experimental demonstra- 
tion we have considered only unitary devices, however, 
it is clear that the method works for non unitary de- 
vices as well. It is also obvious that the method can be 

U M>=n, 8=0.294*] X U [+=0.45lt, 8=-0. 139JE] 

real part 
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FIG. 4: The same experimental characterization as in Fig. |3| 
here for a device made of a series of two optical wave-plates: 
a wave-plate with f> = 0.457T and 9 = — 0.1387T followed by 
another wave-plate with (/> = n and 9 = +0.29-7T. 



used to characterize n-qubit devices — e. g. a controlled- 
NOT gate — in which case we just need multiply by n the 
whole setup, by providing an input entangled state and 
two Pauli detectors for each qubit of the device, with 
the full quantum characterization of the the device ob- 
tained by a joint tomography on all output entangled 
pairs. It is clear that the precision of the method will 
not depend on the particular tested device — whether it 
is unitary or not — and will also be independent on the 
number of qubits. What makes the method particularly 
reliable in the present single-photon polarization encod- 
ing is the fact that all measurements are performed in 
coincidence, making the effect of nonunit quantum effi- 
ciency of detectors negligible, and effectively purifying 
the input entangled state. In a different context — e. g. 
for continuous variables, such as homodyne tomography 
of twin-beams Jt| — quantum efficiency and entanglement 
purity will actually affect the final result: however, the 
quantum tomographic reconstruction can handle all these 
kinds of detection noises below some thresholds |I| , and 
a mixed input state in place of |\&)) work well (but less 
efficiently) as long as the state is faithful, namely it is 
related to a maximally entangled one by an invertible 
map Unfortunately, for the twin-beam homodyne 
tomography [Q, faithfulness requires the knowledge of 
the phase of the pump relative to the local oscillator — a 
feasible but difficult experimental task — whereas in the 
present experiment the form of the entangled input state 
is completely under control, being determined only by 
the orientation of the nonlinear crystal with respect to 
the pump wave-vector and polarization. 



In conclusion, we have given the first demonstration 
of a new tomographic method which allows us to per- 
form a complete characterization of any quantum device, 
exploiting the intrinsic parallelism of quantum entan- 
glement, with a single entangled state playing the role 
of all possible input states in quantum parallel. The 
method works for any generally non unitary multi-qubit 
device, and is particularly reliable in the present context 
of single-photon polarization encoding of the qubit. 
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